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Abstract 

The cross sections for inelastic photoproduction of J/ip and rjj' mesons have 
been measured in ep collisions with the ZEUS detector at HERA, using an inte¬ 
grated luminosity of 38.0 pb“^. The events were required to have 0.1 < z < 0.9 
and 50 < fU < 180 GeV, where z is the fraction of the incident photon energy 
carried by the J/'ip m the proton rest frame and W is the photon-proton centre- 
of-mass energy. The ip' to J/xp cross-section ratio was measured in the range 
0.55 < z < 0.9. The J/b data, for various ranges of transverse momentum, are 
compared to theoretical models incorporating colour-singlet and colour-octet ma¬ 
trix elements. Predictions of a next-to-leading-order colour-singlet model give a 
good description of the data, although there is a large normalisation uncertainty. 
The J/b helicity distribution for z > 0.4 is compared to leading-order QCD 
predictions. 
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1 Introduction 


In the HERA photoproduction regime, where the virtuality of the exchanged photon is 
small, the production of inelastic ^|J mesons, where ifj can be either a J/'0 or a "0', arises 
mostly from direct and resolved photon interactions. In leading-order (LO) Quantum 
Chromodynamics (QCD), the two processes can be distinguished: in direct photon pro¬ 
cesses, the photon couples directly to a parton in the proton; in resolved photon processes, 
the photon acts as a source of partons, one of which participates in the hard interaction. 
Diffractive production, yp —> 'ipN, where is a proton-dissociative state, contributes 
significantly to the inelastic production of ^|J mesons by the direct photon process. 

Direct and resolved photon cross sections can be calculated using perturbative QCD 
(pQCD) in the colour-singlet (CS) and colour-octet (CO) frameworks [1,2,3,4,5,6, 7,8]. 
In the CS model, the colourless cc pair produced by the hard subprocess is identified with 
the physical ^ state. In the CO model, the cc pair emerges from the hard process with 
quantum numbers different from those of the and evolves into the physical 'ijj state by 
emitting one or more soft gluons. At LO, only the photon-gluon-fusion diagram, yp —^ ipg, 
contributes to the direct photon cross section, as shown in Fig. |I|a). Figure |l|b) shows 
the LO diagram for the resolved photon process in the CS framework. A diagram for the 
direct photon process in the CO framework is shown in Fig. |I|c). Although a full next-to- 
leading-order (NLO) calculation of ^ photoproduction is not available for all processes, 
the NLO corrections to the direct photon process, carried out in the framework of the CS 
model, have been calculated [2]. 

The production of ijj mesons has been measured in pp collisions by the CDF collabora¬ 
tion [9,10]. Predictions of the CS model, which for pp collisions exist only at LO in QCD, 
underestimate the data by factors of between 10 and 80. However, after adjustment of 
the corresponding matrix elements, this difference can be accounted for by the CO contri¬ 
butions [4]. Currently, the matrix elements governing the strength of this process cannot 
be calculated, but have to be determined from experiment. Since they are expected to be 
universal, the comparison of the values extracted from ip cross-section measurements in 
different environments constitutes a stringent test of this approach. 

The J/ijj helicity distributions predicted by the CS and CO models have a different de¬ 
pendence on the pt of the J/ip. Furthermore, the dependence of the J/ip polarisation 
on its transverse momentum is sensitive to the virtuality of the initial gluon in the pho¬ 
ton [8]. Results from the CDF collaboration [11] show some discrepancies between the 
helicity measurements and predictions [7] using CO matrix elements extracted from the 
CDF cross-section data. 

The various ip photoproduction processes can be distinguished using the inelasticity vari- 
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able, defined as: 


^ P -v^ 

P ■ q 


( 1 ) 


where P, p^p and q are the fonr-niomenta of the incoming proton, the ip meson and the 
exchanged photon. In the proton rest frame, 2 ; is the fraction of the photon energy 
carried by the ip. Previous HERA data [12,13,14] have shown that the diffractive process 
populates the high-z region, z > 0.9. The direct and resolved photon processes are 
expected to dominate in the regions 0.2 ^ z < 0.9 and z < 0.2, respectively [3]. 

In this study, ip mesons were identified using the decay mode ip —> pPpi~ and were 
measured in the range 50 < IE < 180 GeV, where W is the yp centre-of-mass energy. 
The differential cross sections are given for z > 0.1 and for different regions of transverse 
momentum, pt, of the J/ip. The J/ip helicity distributions in the ranges 0.4 < z < 1 
and 0.4 < 2 ; < 0.9 are presented and compared to model predictions with or without CO 
contributions. 


2 Theoretical models 

2.1 Leading—order Colour Singlet calculations 

The LO matrix element for the photon-gluon-fusion process, as computed in the CS 
framework, is singular for z = 1 and px = 0 [1]. Therefore, the comparisons with these 
theoretical calculations are restricted to the region > 1 GeV. 

Calculations of direct processes at LO in the CS model have been available for some 
time [1]. In this paper, the data are compared to the LO prediction from Kramer et 
ah [2,3] (KZSZ (LO, CS)), including both direct and resolved processes. This calculation 
used the parton density functions (PDFs) GRV94 LO [15] for the proton and GRV LO [16] 
for the photon, the QCD scale parameter, Aq^j^, was set to 0.2 GeV and the factorisation 
and renormalisation scales were set to p = 2mc, where rric, the charm-quark mass, was 
set to 1.5 GeV. Recently, the calculation has been extended to include predictions of the 
J/ip helicity-angle distributions [7] (BKV (LO, CS)). 

In the CS framework the distributions of the J/ip helicity angle have been calculated 
by Baranov [8] for the direct photon process. This calculation uses the /cy-unintegrated 
gluon densities satisfying the BFKL [17] evolution equations. Compared to traditional 
(collinear) parton models, gluons have a transverse-momentum component (or virtuality), 
which results in an increase in the fraction of longitudinally polarised J/ip mesons as pt 
increases. 
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2.2 Next—to—leading—order Colour Singlet calculation 

The NLO corrections to the direct photon process in the CS framework have been calcu¬ 
lated by Kramer et ah [2] (KZSZ (NLO, CS)). This is the only NLO calculation currently 
available for any J/ip production process. The uncertainties in the cross sections arise from 
uncertainties in the non-perturbative QCD parameters. Upper bounds were obtained by 
setting rric = 1.3 GeV and the strong coupling constant, as{Mz), to 0.121 in accordance 
with the MRSTOl {ag t) [18] set of proton PDFs. The lower bounds were obtained by 
setting rric = 1.6 GeV and as{Mz) = 0.117, in accordance with the MRSTOl J.) set of 
PDFs. The dominant uncertainty is due to the variation of the charm-quark mass. For the 
calculation of the cross-section da/dp^, the factorisation and renormalisation scales were 
set to the larger of rric/V^ and -|- /2. For the prediction of the cross-sections 

da/dz, the factorisation and renormalisation scales were set to rric/V2 [2]. 

2.3 Non—relativistic Quantum Chromodynamics calculations 

The LO calculation from Kramer et ah has also been extended to include the GO contri¬ 
butions, (KZSZ (LO, GS-I-GO)), from both direct and resolved photon processes [3]. The 
GO matrix elements were extracted by htting the cross-section da/dp^ for prompt J/ijj 
production measured by GDF [10]. The matrix elements for the hard subprocesses were 
computed at LO, while the GO matrix elements were corrected for initial- and final-state 
gluon radiation by a Monte Garlo (MG) technique [4]. The spread in the predictions is 
due to theoretical uncertainties in the extraction of the GO matrix elements obtained by 
comparing the values extracted by different groups; this spread is often larger than the 
error quoted by each individual group. This calculation has also been extended to predict 
the J/'i/ helicity-angle distributions [7] (BKV (LO, GS-I-GO)). 

A LO calculation by Beneke, Schuler and Wolf [5] (BSW (LO, GS-I-GO)) includes only the 
direct photon process for the GS and GO contributions. Here, the GO matrix elements 
were extracted from measurements by the GLEO collaboration [19] on B meson decays 
to J/'ijj mesons. The matrix elements extracted using the data from GLEO and GDF are 
consistent [3]. This calculation models the so-called shape functions that resum an infinite 
class of GO contributions that are important at high 2 ;. These functions are responsible 
for the decrease of the GO contributions towards z = 1, due to the lack of phase space 
for gluon radiation. This treatment introduces an additional parameter into the model 
called the shape-function parameter which was varied in the range 300 — 500 MeV, based 
on an evaluation [5] of the GLEO data. 

Kniehl and Kramer [6] (KK (LO, GS-I-GO)), like Kramer et ah, have calculated GS and 
GO terms in LO for both direct and resolved photon processes. The GO matrix elements 
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were similarly extracted by fitting the da/dp^ differential cross section for prompt J/'0 
production measured by CDF [10]. The spread in the predictions is due to theoretical 
uncertainties in the extraction of the CO matrix elements. The calculation approximately 
takes into account dominant higher-order (HO) QCD effects and was performed in the 
MS renormalisation and factorisation scheme, using CTEQ4LO [20] and GRV LO as 
the proton and photon PDFs, respectively; the QCD scale parameter was set to 

296 MeV; common factorisation and renormalisation scales were used and were set to 
/i = -|- py with 2mc = rrij/jp. 

3 Experimental conditions 

The data were collected during the 1996 and 1997 running periods, when HERA operated 
with protons of energy Ep = 820 GeV and positrons of energy Eg = 27.5 GeV, and 
correspond to an integrated luminosity of 38.0 ± 0.6 pb~^. This represents more than a 
tenfold increase with respect to the previous ZEUS analysis [12]. A detailed description 
of the ZEUS detector can be found elsewhere [21,22]. A brief outline of the components 
that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [23], which operates 
in a magnetic held of 1.43 T provided by a thin superconducting coil. The CTD consists 
of 72 cylindrical drift chamber layers, organized in 9 superlayers covering the polar an¬ 
gle^ region 15° < 9 < 164°. The transverse-momentum resolution for full-length tracks is 
= 0.0058pT © 0.0065 © 0.0014/^^, with pt in GeV. 

The high-resolution uranium-scintillator calorimeter (GAL) [24] consists of three parts: 
the forward (FGAL), the barrel (BGAL) and the rear (RGAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec¬ 
tion (EMG) and either one (in RGAL) or two (in BGAL and FGAL) hadronic sections 
(HAG). The smallest subdivision of the calorimeter is called a cell. The GAL energy res¬ 
olutions, as measured under test beam conditions, are a{E)/E = 0.18/y/E for electrons 
and a{E)/E = 0.35/y/E for hadrons {E in GeV). The timing resolution of the GAL is 
better than 1 ns for energy deposits greater than 4.5 GeV. 

The muon system consists of tracking detectors (forward, barrel and rear muon chambers: 
FMUON [21], B/RMUON [25]), which are placed inside and outside a magnetized iron 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the “forward direction”, and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as ?7 = — In (tan |), where the polar angle, 0, is measured with respect to the proton beam 
direction. 
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yoke surrounding the CAL and cover polar angles from 10° to 171°. The barrel and 
rear inner muon chambers cover polar angles from 34° to 135° and from 135° to 171°, 
respectively. 

The luminosity was determined from the rate of the bremsstrahlung process e~^p —>■ 
where the photon was measured by a lead-scintillator calorimeter [26] located at Z = 
— 107 m. 


4 Event selection and ijj reconstruction 

The -0 —>■ candidates were selected using a three-level trigger system [21]. In the 

hrst-level trigger (FLT), the barrel and rear inner muon chambers, BMUI and RMUI, 
were used to tag the muons from -0 decays by matching segments in the muon chambers 
with tracks in the CTD, as well as with energy deposits in the CAL consistent with the 
passage of a minimum ionising particle (m.i.p.). Events satisfying this regional matching 
and having tracks in the CTD pointing to the nominal interaction vertex were selected. 

In the second-level trigger (SLT), the total energy in the calorimeter {E = SjEj) and the 
Z component of the momentum [pz = 'ZiEi cos 6i) were calculated. The sums run over all 
calorimeter cells i with an energy, Ei, and polar angle, 6i, measured with respect to the 
nominal vertex. To remove proton-gas interactions, events with the ratio pzjE greater 
than 0.96 were rejected. The cosmic-ray background was partially rejected at the SLT 
by using the time differences of energy deposits in the upper and the lower halves of the 
BCAL. 

In the third-level trigger (TLT), a muon candidate was selected when a track found in 
the CTD matched both a m.i.p. in the CAL and a track in the inner muon chambers. 
An event containing a muon candidate in the rear (barrel) region was accepted if the 
momentum (transverse momentum) of the CTD track exceeded 1 GeV. 

In the offline analysis, the TLT algorithm was again applied to the results of the full 
event reconstruction. In addition, the tracks corresponding to the two muons from the 
^|J decay had to satisfy several criteria. One track was matched to both a m.i.p. cluster 
in the CAL and a track in the inner muon chambers. This track was required to have a 
momentum greater than 1.8 GeV if it was in the rear region, or a transverse momentum 
greater than 1.4 GeV if in the barrel region. The other muon track was matched to a 
m.i.p. cluster in the GAL and was required to have a transverse momentum of greater 
than 0.9 GeV. Both tracks were restricted to the pseudorapidity region |? 7 | < 1.75. To 
reject cosmic rays, events in which the angle between the two muon tracks was larger 
than 174° were removed. Events were also required to have an energy deposit greater 
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than 1 GeV in a cone of 35° around the forward direction (excluding calorimeter deposits 
due to the decay muons). Elastically produced -0 mesons were thus excluded. 

The kinematic region considered was defined by the inelasticity variable given in Eq. (^), 
and by the photon-proton centre-of-mass energy: 

fE2 = (P + g)2; 

W and 2 ; were computed from 

= ‘2E^{E-vz) 

and 

_ (P - Vz)'i> 

E-vz ’ 

where E - pz = {E - pz)had + {E - pz)^p- The quantity {E - pz)ha.d is the sum over 
the hadronic hnal state, calculated using all CAL cells excluding those belonging to the 
decay-muon clusters; {E — pz)^) was calculated using the -0 decay tracks measured by the 
CTD. 

The events were required to have E — pz < 20 GeV, which restricts W to be less than 
180 GeV and < 1 GeV^, with a median value of ~ 10“^ GeV^. The elimination of 
deep inelastic scattering events was independently conhrmed by searching for scattered 
positrons in the CAL [27]; none was found. As the analysis uses only the B/RMUON, 
the polar angle coverage of these detectors restricts W to be greater than 50 GeV. 


5 Monte Carlo models 

The production of ijj mesons from direct interactions was simulated using the HER- 
WIG 5.8 [28] MG generator, which generates events according to the LO diagrams of 
the photon-gluon-fusion process, jg —*• 'ijjg, as computed in the framework of the CS 
model. The hadronisation process is simulated by the cluster model [29]. Events were 
generated in the range of starting from the kinematic limit (~ 10“^° GeV^) up to 
10 GeV^. Events were generated for 2 < 0.95 to avoid a singular phase-space region. The 
GRV94 LO PDF for the proton was used. The HERWIG MG sample was reweighted in 
Pt and W to the data. 

The production of ^ mesons from resolved photon interactions was simulated using the 
PYTHIA 6.146 [30] MG generator (resolved photon interactions for ^|J production are not 
implemented in HERWIG). The GRV94 LO and GRV LO PDFs were used for the proton 
and photon, respectively. The matrix elements for the resolved photon processes were 
computed at LO in the framework of the CS model. 
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Diffractive production of V' mesons with proton dissociation was simulated with the EP- 
SOFT [31] MC generator, which has been tuned to describe such processes at HERA [32], 
Data in the region 0.9 < z < 1 were used to determine the dependence of the cross sec¬ 
tion on the invariant mass of the dissociative system, on the photon-proton centre-of-mass 
energy and on the of the J/ip meson. 


6 Signal determination and cross-section calculation 

The invariant-mass spectra of the muon pairs measured for 50 < PE < 180 GeV and three 
ranges are shown in Fig. The J/ip is clearly seen in all ranges. The higher ranges 
(Figs. Hb) and |]c)) also show a ip' peak. 

To estimate the number of events in the signal regions (the mass ranges 2.9 to 3.25 GeV 
and 3.6 to 3.8 GeV for the J/ip and ip', respectively), an accurate description of the 
combinatorial background was necessary. This was estimated by htting the invariant- 
mass distribution of the data outside the corresponding windows of the ip mesons, using 
a function which is the product of a second-order polynomial and an exponential. 

The data were corrected bin by bin for geometric acceptance, detector, trigger and re¬ 
construction efficiencies, as well as for detector resolution. The correction factor, as a 
function of an observable (P in a given bin i, is Ci{0) = Af®'^(C>)/A')'®'^((P). The variable 
^gen(^) .g number of generated events and is the number of reconstructed 

events passing the selection requirements detailed in Section ^ Both numbers were com¬ 
puted using the MG generators described in Section In this W range and the three 
regions of z, 0.1 < z < 0.4, 0.4 < .2 < 0.9 and 0.9 < .2 <1, the acceptance (dehned as 
1/Ci{0)) was typically 30% and always above 10%. 

For 0.9 < z < \, the events are largely diffractive. Therefore, the analysis of inelastic 
production was restricted to the region 0.1 < z < 0.9. The remaining contamination 
from diffractive processes was estimated by htting the relative fractions of HERWIG and 
EPSOFT MG event samples to the data. A £t was performed to the inelasticity 
distribution in the region 0.4 < z < 1 and three Pt ranges: 0<p7’<l, 1 < Pt < ‘2 and 
Pt > 2 GeV. Figure ^shows, for events in the region 50 < PE < 180 GeV and 0.4 < z < 1, 
that the resulting mixture of 56% HERWIG and 44% EPSOFT gives a reasonable descrip¬ 
tion of the relevant J/ip event observables. For 0.4 < 2 ; < 0.9, with no pr cut, the diffrac¬ 
tive contribution, as estimated with the EPSOFT MG, was 17%, concentrated at low px- 
The diffractive contribution was subtracted bin by bin for all cross-section measurements. 

Resolved photon processes are also present in the region of the cross-section measurement. 
In the region 0.1 < z < 0.9, the resolved photon component was estimated by htting the 
relative fractions of direct and resolved photon events in the MG samples to the inelasticity 


7 


distribution in the data. The fraction of resolved photon events is 5%, reaching up to 
50% for 0.1 < 2 ; < 0.4, in good agreement with theoretical expectations [3]. In this low 2 ; 
range the contribution from B decays can be as large as 25%; it is negligible at higher z. 
The photoproduction cross section was obtained from the measured electron-proton cross 
section by dividing by the integrated effective photon flux [12]. 


7 Systematic uncertainties 

A detailed study of possible sources of systematic uncertainties was carried out for all 
measured differential cross sections. The following sources were considered: 

• muon-chamber efficiencies: the BMUI and RMUI muon-chamber efficiencies were ex¬ 
tracted from the data using muon pairs coming from elastic J/ip events and from the 
process 77 —> ii~ ■ For > 1-4 GeV, the product of the geometrical and chamber 
efficiency for the BMUI is greater than 30%, reaching 60% at high transverse momen¬ 
tum. For p > 1.8 GeV, this product for the RMUI is greater than 45%, reaching 70% 
at high momentum. The associated uncertainty of about ±7% is independent of the 
phase-space region; 

• analysis cuts: this class comprises the systematic uncertainties due to the uncertainties 
in the measurement of momentum, transverse momentum and pseudorapidity of the 
muon decay tracks. Each cut was varied within a range determined by the resolution in 
the appropriate variable. The pseudorapidity contribution gave an uncertainty below 
± 1 %, while variation of the appropriate momentum or transverse momentum cut gave 
a ±1.5% contribution; 

• GAL energy scale: GAL energy measurements were used in the W and ^ reconstruc¬ 
tion. This leads to a systematic uncertainty in the measured cross section due to 
the ±3% uncertainty on the energy scale of the GAL. This effect was investigated by 
varying the quantity {E — pz)ha.d by ±3% in the MG sample, leading to a variation of 
the cross sections of ± 6 % for 0.1 < z < 0.4. Integrated over z, the effect is typically 
below ± 2 %; 

• GAL energy resolution: the W and z resolutions are dominated by the GAL energy 
resolution through the quantity {E — pz)had- The {E — pz)ha.d resolution in the MG 
was smeared event by event by ±20%. This estimated any possible mismatch between 
the {E — pz)hsid resolution in the data and MG simulation, giving an uncertainty of 
±3%; 

• diffractive subtraction: the fraction of HERWIG and EPSOFT MG events, hxed by 
the htting procedure described in Section ^ is known to a precision limited by the 


number oi J/^p events in the data and the process modelling. All possible fractions 
giving a iii the interval [Xmin; Xmin + t] were considered and the largest change in 
the cross section was quoted as the systematic uncertainty. This gave an uncertainty 
which was at most ±2% at high z and low where the diffractive contribution peaks; 

• diffractive simulation: the EPSOFT MC simulation parameters were varied within 
ranges allowed by the comparison between the data and the EPSOFT MC simulation 
in the region 0.9 < z < 1. The fraction of HERWIG and EPSOFT MC events was 
re-evaluated. This gave an uncertainty which was at most ±2.5% at high 2 ; and low 
Pt] 

• Pt and W spectra: the px and W spectra of the J/ip meson in the HERWIG MC 
simulation were varied within ranges allowed by the comparison between the data and 
the simulation and the correction factors re-evaluated. This gave an uncertainty of 
± 2 %; 

• helicity parameterisation uncertainty: in the HERWIG MC, the helicity parameter a 
is set to 0. As a systematic check, the helicity in the HERWIG MC was reweighted 
according to the upper and lower limits of error in the measured distribution and the 
correction factors re-evaluated. This gave an uncertainty of ±5%, independent of the 
phase-space region. 

All of the above individual sources of systematic uncertainty were added in quadrature. 
The following sources resulted in an overall shift of the cross section and were therefore 
treated separately: 

• the integrated luminosity determination gave an uncertainty of ±1.6%; 

• the branching ratio of J/ip —> gave an uncertainty of ±1.7% [33]. 


8 Results 

8.1 Total cross-section measurement at high ^ 

A cross-section measurement in the region 2 ; > 0.9 is particularly interesting because the 
CO mechanism is expected to contribute signihcantly at high 2 ; [3], whereas for 2 ; < 0.9 the 
sensitivity to this production mechanism is reduced. In particular, at large . 2 , the impact 
of the soft-gluon emission on the hadronisation of CO cc pairs is not well understood [3, 
5,34]. In the region 0.9 < z < 1, the 2 ; resolution is comparable to the width of the 2 ; 
interval and the diffractive process is dominant; hence the separation of the direct and 
diffractive components is not reliable. Therefore, only the visible cross section in the 
region 0.9 < z < 1, including the diffractive component, is given. Furthermore, due to 
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the requirement of an energy deposit in the direction of the outgoing proton, necessary to 
remove the elastic component, only diffractive events with a high-enough invariant mass, 
Mtv, of the hnal-state hadronic system were included here. Monte Carlo studies showed 
that the requirement of an energy deposit exceeding 1 GeV in a 35° cone around the 
outgoing proton direction corresponds to a threshold in M^v of 4.4 GeV, above which all 
correction factors were independent of M^- The following cross sections, corresponding 
to the phase-space region dehned by 50 < hh < 180 GeV, 0.9 < z < 1, > 4.4 GeV 

were obtained: 


(^j/tpiPT > 0 GeV) 
(^j/ip{PT > 1 GeV) 
^j/pjiPT > 2 GeV) 


45.7 ± 1.3 (stat.) (syst.) nb; 

24.5 ± 0.9 (stat.) (syst.) nb; 

6.5 ± 0.5 (stat.) Ig;® (syst.) nb. 


The uncertainties on the value of the Mm threshold and the GAL energy resolutions are 
of similar importance and dominate the systematic uncertainty. 


8.2 Measurement of production 


The production of xjj' with subsequent decay to J/'0 has been measured using the rates of 
■0' —»■ pi~ and J/ip ^ ■ The ip' to J/pj cross-section ratio was determined in the 

region 0.55 < z < 0.9 with no px cut on the ip mesons. The range 0.4 < z < 0.55 was not 
included, because it has a large non-resonant background. The xp' to J/pj cross-section 
ratio was computed in bins of W and 2 ; from 

aiPP') _Nf Cf BM r Nf Cf BM 
CT*(J/V') “ ^ ' V ^ ^ ’ 

where, for the considered bin i, is the number oi J/ip pip') events observed, 

C}^ is the correction factor (see Section ||) computed using HERWIG MG J/pj {pj') 

events, Br^ (Br^') is the J/pj pip') muonic branching ratio and Br' is the pj' ^ J/p} X 
branching ratio. The values used were BM = (5.88 ± 0.10)%, Br^' = (0.70 ± 0.09)% and 
Br' = (55.7 ±2.6)% [33]. 

With this technique, the cross-section ratio is corrected for the pj' —> J/pj {—>■ p~) X 

cascade decay. The results are shown in Fig. | and listed in Table |l|; all cross-section 
ratios are consistent with being independent of the kinematic variable, as expected if 
the underlying production mechanisms for the J/pj and pj' are the same. For the range 
0.55 < 2 ; < 0.9 and 50 < IV < 180 GeV, the pj' to J/pj cross-section ratio is 


(T{J/pj) 


0.33 ± 0.10 (stat.)lQ;Q2 (syst-)) 


( 2 ) 


10 










in agreement with the expectation of the LO CS model of 0.24 [2]. 

Even though the NLO corrections to the CS model for J/'ij) production are known to be 
large, similar large NLO corrections are expected to affect the rate of production [35]. 
Hence the "0^ to J/V' cross-section ratio at NLO is not expected to differ signihcantly from 
that at LO. From the cross-section ratio and the iIj' —*• J/'0X branching ratio [33], it is es¬ 
timated that the observed cross section for mesons is increased by (18.4±5.6 (stat.))% 
due to J/'ll! mesons originating from '0' cascade decays. This is consistent with the ex¬ 
pected value of 15% [2], which has been added to all predictions of the differential 
cross sections presented in this paper. 


8.3 Measurement of inelastic J/t/; cross sections 


The differential cross-sections da/dz for three different regions in pt and 0.1 < z < 0.9 
are shown in Figs. BS and 1^ and listed in Tables and |^. All data sets show a cross 
section increasing with 2 ;. The differential cross-section da/dp^, measured in the region 
0.4 < < 0.9, is shown in Fig. and listed in Table |^. The measurement extends to 

p"^ ~ 24 GeV^, where the cross section has fallen by two orders of magnitude. The cross- 
section has a function of W, for > 1 GeV and 0.4 < 2 ; < 0.9, is given in Table |] and 
shown in Fig. |^. The differential cross-section da/dy in the region 0.4 < z < 0.9 and 
Pt > 1 GeV is shown in Fig. 0 and listed in Table |^. The rapidity, p, of the J/V' is given 
by 



E + pz 
E-pz' 


where E and pz are the energy and Z component of the momentum of the J/'0, respec¬ 
tively. 


The helicity distribution can be parameterised as 


dN 

d cos 6* 


oc 1 -|- a cos^ 9*. 


(3) 


The helicity-parameter a was determined by reweighting the HERWIG MG dN/d cos 9* 
generator-level distribution according to Eq. for different values of a. The for fho 
dN/d\ cos 0*1 distribution in data and MG was then calculated for each value of a in the 
MG and the minimum gave the central value of a. The procedure was repeated for 
each Pt bin in the range 1 < pt < 5 GeV. The systematic uncertainties were negligible 
with respect to the error obtained from the ht. 

Figure ^ shows the measured parameter a plotted as a function of the pt of the J/V'- 
In Fig. |IT|a) and b), the quantisation axis is chosen to be the opposite of the incoming 
proton direction in the J/'ijj rest frame, 9* is the opening angle between the quantisation 
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axis and the /i’*' direction of flight in the J/ip rest frame and a is the helicity parameter. 
This frame is known as the “target frame”. The parameter a was determined in bins of 
Ptj for pt > I GeV and 0.4 < z < 1 (Fig. pT|a)). The parameter a was also measured in 
the range OA < z < 0.9 (Fig. 13))). where the diffractive contamination is reduced. The 
values of the helicity-parameter a are summarized in Table p. The CDF helicity-angle 
dehnition was also used, where the quantisation axis was dehned as the J/ijj direction of 
flight in the ZEUS coordinate system; this frame is known as the “helicity basis” [11]. The 
parameter a measured in this frame is shown in Fig. pT|c) and d) and listed in Table |^. The 
values a = —1 and a = +1 correspond to fully longitudinal and transverse polarisation, 
respectively. Within the large experimental uncertainties, the data are consistent with a 
trend from transverse to longitudinal polarisation with increasing px. 

These results are consistent with the previous ZEUS measurements [12], but have an 
improved precision and extend to higher p\. Recently, the HI collaboration has also 
published a measurement of inelastic J/V' production [14] showing similar features to 
those presented here. 


8.3.1 Comparison with leading—order Colour Singlet calcula¬ 
tions 

Figure shows a comparison of the KZSZ (LO, CS) predictions with the data. For 
~ 1 GeV^, the prediction underestimates the data by a factor of about two, although 
this is within the range of the theoretical uncertainties. For higher px values, the cal¬ 
culation falls increasingly below the data. At p^ ~ 20 GeV^, the LO GS prediction 
undershoots the data by a factor of ~ 20. 

The prediction of BKV (LO, GS) for the helicity-parameter a as a function of pt in the 
target frame is shown in Figs. pT|a) and b). The prediction lies somewhat below the data 
at low Pt and somewhat above at high pt, although the data have large statistical uncer¬ 
tainties. This general trend appears for both 0.4 < z < 1 and 0.4 < z < 0.9 regions. The 
prediction from Baranov also lies below the data at low px, but gives a good description 
of the data for > 1.6 GeV. 

The LO GS prediction from Baranov is also shown in Fig. pT|c) and d) compared to the 
data for the helicity-base frame. The GRV prediction was obtained by adding [36] a kx 
dependence to the GRV (collinear) gluon density. The KMS prediction was obtained from 
the /cT-unintegrated gluon density [37]. The data is reasonably well described by both 
predictions. 
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8.3.2 Comparison with the next—to—leading—order Colour Sin¬ 
glet calculation 

The KZSZ (NLO, CS) prediction is compared to the data for da/dp^ in Fig. ||. The 
prediction is not reliable in the pt ^ 0 limit [2] and hence is only shown for pt > 1 GeV. 
The predicted shape, controlled by QCD glnon radiation, is in good agreement with the 
data. The normalisation of the predicted cross sections is sensitive to the assnmed valnes 
of the non-pertnrbative QCD input parameters, such as the mass of the charm quark and 
the value of Aqcd- The uncertainty resulting from the variation of the renormalisation and 
factorisation scales is small by comparison [35]. The overall normalisation of the data is 
well described by the calculation, although the prediction suffers from large uncertainties. 
The discrepancy between the data and the LO CS prediction can therefore be explained 
by large NLO corrections. 

The prediction is also compared to the cross-section as a function of W in Fig. |^. The 
prediction again suffers from large theoretical uncertainties but describes the shape and 
normalisation of the data. Similar conclusions can be drawn from the comparison with 
the differential cross-sections da/dz shown in Figs. 

8.3.3 Comparison with non—relativistic Quantum Chromody¬ 
namics calculations 

The inelasticity distributions in Figs. are compared with different predictions includ¬ 
ing CO matrix elements extracted by htting independent data sets. In Fig. H, da/dz 
for pt > 1 GeV is compared with the KZSZ (LO, CS-I-CO) and KK (LO, CS-I-CO) cal¬ 
culations. The rise of the predicted cross section for z < 0.1 is due to resolved photon 
processes. Within the large theoretical uncertainties, the prediction KZSZ (LO, CS-I-CO) 
gives a good description of the data. The KK (LO, CS-I-CO) result lies signihcantly below 
the data, but describes the shape reasonably well. 

In Fig. p, the differential cross-section da/dz for px > 2 GeV is compared with the 
BSW (LO, CS) and BSW (LO, CS-I-CO) calculations [5]. The CS prediction clearly lies 
below the data, whilst the inclusion of the CO terms gives a better description. The spread 
in the prediction, which is largest at high z, is due to the uncertainty on the value of the 
shape-function parameter. The overall shape of the spectrum is weakly dependent on the 
CO matrix elements, which primarily affect the global normalisation of the spectrum. 

The KK (LO, CS-I-CO) prediction is compared to the cross-section da/dy in Fig. ^ for 
the kinematic region 50 < hF < 180 GeV, 0.4 < z < 0.9 and pr > 1 GeV. The calculation 
includes both direct and resolved photon processes, but the resolved photon contribution 
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is negligible in the selected phase-space region (dne to the lower 2 ; cnt). The predicted 
cross section falls well below the data, although the shape is reasonably reproduced, as 
shown when the prediction is multiplied by a factor of three. 

In Figs. pT]a) and b), the BKV (LO, CS-I-CO) prediction is compared to the data for the 
helicity-parameter distribution as a function of pt- The spread in the prediction is due 
to theoretical uncertainties on the values of the CO matrix elements. In the currently 
accessible px range, the CS plus CO predictions are similar to those of the CS only, 
although the prediction from the CS model rises with p-r, while the CS plus CO prediction 
decreases slightly. 


9 Conclusions 

Cross sections of inelastic J/'0 photoproduction have been measured and compared with 
LO and NLO QCD predictions. The LO CS prediction does not describe the spectrum. 
A NLO QCD calculation in the framework of the CS model, including only the direct 
photon process, gives a good description of the p^ and 2 ; differential cross sections and of 
the cross section as a function of W. However, given the large theoretical uncertainties 
affecting the NLO calculation, it is currently not possible to constrain the size of the CO 
contributions. Furthermore, LO calculations including CO contributions, as determined 
from pp data, describe the data, albeit with large theoretical uncertainties. Although the 
helicity distribution at high p^ is sensitive to the underlying production mechanism, the 
data are unable to distinguish between the two mechanisms. These results agree with 
the measurements recently published by the HI collaboration after taking into account, 
by MC extrapolation, the small (< 10 %) normalisation differences due to the different 
phase-space regions probed by the two experiments. 
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Pt 

(GeV) 

<Pt> 

(GeV) 

a{4')/a{J/4) 

0.0 - 1.0 

0.5 

0.17±0.12 

1.0 - 1.75 

1.38 

0.35±0.17 

1.75 - 5.0 

3.38 

0.26T0.15 

IT 

(GeV) 

< IT > 

(GeV) 

a{4')/a{J/4) 

50. - 85. 

67.5 

0.28±0.10 

85. - 110. 

97.5 

0.26±0.14 

no. - 180. 

145. 

0.32T0.20 

z 

< z > 

a{4')/a{J/4) 

0.55 - 0.7 

0.625 

0.39±0.26 

0.7 - 0.8 

0.75 

0.19±0.12 

0.8 - 0.9 

0.85 

0.30±0.10 


Table 1: Cross-section ratios between 4’' and J/4 as function of pt, W and z 
variables. These ratios are measured in the kinematical region 50 < IT < 180 GeV 
and 0.55 < < 0.9. The uncertainties are statistical only. 


z 

< 2 ; > 

bkg 

da/dz 



% 

(nb) 

0.10-0.40 

0.28 

0 . 

2Q.Q±V2tll 

0.40-0.55 

0.47 

0 . 

45.3±5.5;i^ 

0.55-0.70 

0.62 

3. 

76.3±4.i;®;^ 

0.70-0.80 

0.75 

9. 

96.7±4.4l^^/ 

0.80-0.90 

0.85 

34. 

97.6±3.5l^^8^ 


Table 2: Differential cross-section da/dz measured for 50 < IT < 180 GeV 
and no pr cut. In the quoted cross sections, the first uncertainty is statistical and 
the second is systematic. Overall normalization uncertatinties due to the luminos¬ 
ity measurement (±.1.Q%) and to the J/fj decay branching ratio (±1.7%) are not 
included in the systematic error. The column labelled bkg gives, in each bin, the 
percentage of diffractive background subtracted from the data. 


17 


































< z > 

bkg 

% 

da/dz 

(nb) 

0.10-0.40 

0.28 

0 . 

16.3±3.4+^;^ 

0.40-0.55 

0.47 

0 . 

27.1±4.4+|° 

0.55-0.70 

0.62 

2 . 

49.0±3.3lt9 

0.70-0.80 

0.75 

8 . 

66.2±3.7+®;® 

0.80-0.90 

0.85 

28. 

68.9±3.1;|3 


Table 3: Differential cross-section da/dz measured for 50 < IT < 180 GeV and 
Pt > I GeV. In the quoted cross sections, the first uncertainty is statistical and 
the second is systematic. Overall normalization uncertatinties due to the luminos¬ 
ity measurement (±1.6%) and to the J/fi decay branching ratio (±1.7%) are not 
included in the systematic error. The column labelled bkg gives, in each bin, the 
percentage of diffractive background subtracted from the data. 


2 : 

< ^ > 

bkg 

da/dz 



% 

(nb) 

0.10-0.40 

0.28 

0 . 

8 .0±2.2+J;° 

0.40-0.55 

0.48 

0 . 

12.1±2.5+?| 

0.55-0.70 

0.62 

0 . 

20.4±2.2+2;4 

0.70-0.80 

0.75 

6 . 

23 . 0 ± 2 . 1+|2 

0.80-0.90 

0.85 

16. 

31.3±2.2;^;^ 


Table 4: Differential cross-section da/dz measured for 50 < IT < 180 GeV and 
Pt > 2 GeV. In the quoted cross sections, the first uncertainty is statistical and 
the second is systematic. Overall normalization uncertatinties due to the luminos¬ 
ity measurement (±1.6%) and to the J/fi decay branching ratio (±1.7%) are not 
included in the systematic error. The column labelled bkg gives, in each bin, the 
percentage of diffractive background subtracted from the data. 
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Pt 

<Pt> 

bkg 

da/dp^ 

(GeV2) 

(GeV^) 

% 

(nb/GeV^) 

0. 

- 1. 

0.46 

22. 

12.3±0.5lg;^ 

1. 

- 2. 

1.45 

18. 

7.5±0.4+°J 

2. 

~ 3. 

2.45 

16. 

4.8±0.3l[ji 

3. 

- 4. 

3.46 

15. 

2.7±0.3l|JJ 

4. 

~ 5. 

4.47 

14. 

2.3±0.3l|J;^ 

5. 

- 6. 

5.40 

10. 

1.73±0.23+°;^® 

6. 

~ 7. 

6.45 

3. 

1.56±0.21+°;i^ 

7. 

~ 8. 

7.55 

0. 

0.94±0.161^6 

8. 

- 9.5 

8.78 

0. 

0.71±0.12l°;?® 

9.5 

- 11.5 

10.48 

0. 

0.50±0.08ino 

11.5 

- 15.5 

13.53 

0. 

0.24±0.04l|]:g^ 

15.5 

- 21. 

18.02 

0. 

0.15±0.02l°;°^ 

21. 

- 30. 

24.22 

0. 

0.043±0.009l°;°J^ 


Table 5: Differential cross-section da/dpf measured for bO < W < 180 GeV and 
0.4 < 2 < 0.9. In the quoted cross sections, the first uncertainty is statistical and 
the second is systematic. Overall normalization uncertatinties due to the luminos¬ 
ity measurement (±1.6%) and to the J/fj decay branching ratio (±1.7%) are not 
included in the systematic error. The column labelled bkg gives, in each bin, the 
percentage of diffractive background subtracted from the data. 


W 

< W > 

bkg 

a 

(GeV) 

(GeV) 


(nb) 

50- 70 

61.3 

27. 

22.0±1.5l|^ 

70- 90 

79.9 

14. 

24.1±1.4l|^ 

90 - 110 

100.1 

10. 

29.1±1.9l2| 

0 

i 

0 

124.8 

9. 

28.9±2.0l|^ 

140 - 180 

157.2 

4. 

29.6±3.2l|® 


Table 6: Cross section versus W measured for 0.4 < 2 ; < 0.9 and pt > 1 GeV. 
In the quoted cross sections, the first uncertainty is statistical and the second is 
systematic. Overall normalization uncertatinties due to the luminosity measure¬ 
ment (±1.6%) and to the J/fi decay branching ratio (±1.7%) are not included in 
the systematic error. The column labelled bkg gives, in each bin, the percentage of 
diffractive background subtracted from the data. 


19 
































y 

<y> 

bkg 

% 

da/dy 

(nb) 

- 1 . 6 -- 1.2 

-1.32 

7. 

3.3±0.5lgi 

- 1 . 2 -- 0.8 

-0.97 

9. 

7.9±0.7+°| 

- 0.8 - - 0.4 

-0.59 

9. 

ll.l±0.8li;^ 

- 0.4 - 0 . 

-0.20 

10. 

10.5±0.7li;2 

0 . - 0.4 

0.18 

14. 

10.6±0.7+};J 

0 . 4 - 0.8 

0.59 

18. 

10.9±0.7ll:? 

0 . 8 - 1.2 

0.92 

17. 

8.1±0.9l};)j 


Table 7: Differential cross-section da/dy measured for 50 < IT < 180 GeV, 
0.4 < z < 0.9 and pt > I GeV. In the quoted cross sections, the first uncertainty 
is statistical and the second is systematic. Overall normalization uncertatinties 
due to the luminosity measurement (±1.6%) and to the J/fi decay branching ratio 
(±1.7%) are not included in the systematic error. The column labelled bkg gives, 
in each bin, the percentage of diffractive background subtracted from the data. 


Pt 

(GeV) 

<Pt> 

(GeV) 

a (0.4 < z < 1) 

a (0.4 < z < 0.9) 

1.0 - 1.2 

1.10 

1 19 + 0.72 

1 1 9+0*93 
-‘-•-‘-^-0.74 

1.2 - 1.4 

1.30 


0.82l°;?3 

1.4 - 1.6 

1.50 

0 76+°-^2 

' 'J-0.58 


1.6 - 1.9 

1.74 

ri qo+0.63 
'^•^'^-0.51 

U-'JU o.63 

1.9 - 2.4 

2.13 

-0 09'''°'^^ 

'^■^■‘--0.44 

2.4 - 3.4 

2.81 

U-UO g 37 

-0 

'^•"^'-0.32 

3.4- 5.0 

4.06 

U-^+_0.48 

-0 

59 


Table 8: J/f> helicity parameter a as a function of pt measured in the target 
frame for 50 < IT < 180 GeV and 0.4 < 2 ; < 1(0.9). The uncertainties are due to 
the total experimental uncertainties. 
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Pt 

(GeV) 

<Pt> 

(GeV) 

a (0.4 < z < 1) 

a (0.4 < 2 ; < 0.9) 

1.0 - 1.2 

1.10 

-0.33 

fi 07 + 0.52 
'J-'J'0.42 

1.2 - 1.4 

1.30 

o.29 

-0 05’''°''^^ 

U-Ud 0.35 

1.4 - 1.6 

1.50 

-0 07 

-0.29 

n e:e:+0.62 

U.00-0,48 

1.6 - 1.9 

1.74 

-0.09 

-0 1 

'^■■‘■'J-0.38 

1.9 - 2.4 

2.13 

-0.28 

-0.491°:^^ 

2.4 - 3.4 

2.81 

fi 97 +0.52 

Q42 

-0.58l°i^ 

3.4- 5.0 

4.06 

0.39 tlit 

0.74ll:i 


Table 9: J/'0 helicity parameter a as a function of pt measured in the helicity 

basis for 50 < IT < 180 GeV and 0.4 < z < 1(0.9). The uncertainties are due to 
the total experimental uncertainties. 
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Figure 1: a) The direct photon process at leading order in the CS framework; b) 
the resolved photon process in the same framework; c) the direct photon process in 
the CO framework. 
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Figure 2: The invariant-mass spectrum measured in the region 

50 <W < 180 GeV for a) 0.1 < z < 0.55, b) 0.55 < z < 0.9 and c) 0.9 
< z < 1. The signal regions are shown as the shaded bands and the background as 
the continuous line. 
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Figure 3: Number of events reconstructed in the kinematic region z > O .4 and 
50 < W < 180 GeVplotted against a) J/fj px, b) W c) inelasticity, z and d) J/if 
polar angle, 6. The data distributions are shown as the points with statistical errors 
only. The simulated EPSOFT diffractive proton-dissociation background is shown 
as the dotted lines. The solid lines show the prediction of the sum of the HERWIG 
and EPSOFT generators. The combined MO has been area normalised to the data. 
The HERWIG MO sample was reweighted in px and W to the data. 
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Figure 4: -0' to J/'0 cross-section ratios as a function of a) px, b) W, and c) 

z, measured for 50 < W < 180 GeV and 0.55 < z < 0.9. The uncertainties are 
statistical only. 
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Figure 5: Differential cross-sections da/dz for 50 < PF < 180 GeV and different 
Pt selections: no pr cut (squares), pt > 1 (circles), and pr > 2 GeV (trian¬ 
gles). The inner error bars show the statistical uncertainty; the outer bars show the 
statistical and systematic uncertainties added in quadrature. 
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Figure 6: J/V' differential cross-section da/dp^ for 50 < W < 180 GeV and 
0.4 < z < 0.9. The inner error bars show the statistical uncertainty; the outer bars 
show the statistical and systematic uncertainties added in quadrature. The data 
points are compared to the prediction KZSZ (NLO, CS) (shaded band) including 
only the direct photon process. The spread in the prediction is due to uncertainties 
on the charm-quark mass and on the QCD scale parameter, Aqcd- The dotted line 
represents the LO prediction KZSZ (LO, CS). A 15 % contribution has been added 
to the predictions to account for J/^f> mesons originating from if' cascade decays. 
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Figure 7: Cross section as a function of W for pr > 1 GeV and O.f < z < 0.9. 
The inner error bars show the statistical uncertainty; the outer bars the statistical 
and systematic uncertainties added in quadrature. The shaded band shows the pre¬ 
diction KZSZ (NLO, CS). The spread in the prediction is due to uncertainties on 
the charm-quark mass and on the QCD scale parameter, Aqcd- A 15 % contribu¬ 
tion has been added to the prediction to account for J/fj mesons originating from 
xjj' cascade decays. 
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Figure 8: Differential cross-section da/dz for 50 < PF < 180 GeV and 

Pt > 1 GeV (points). The inner error bars show the statistical uncertainty; the 
outer bars show the statistical and systematic uncertainties added in quadrature. 
The shaded band shows the prediction KZSZ (NLO, GS), including only the direct 
photon process. The spread in the prediction is due to uncertainties on the charm- 
quark mass and on the QGD scale parameter, Kqcd- The solid lines show the 
prediction of the KZSZ (LO, GS-hGO) calculation performed including both direct 
and resolved photon processes. The spread in the predictions is due to theoretical 
uncertainties in the extraction of the GO matrix elements. The KK (LO, OS-hOO) 
prediction is also shown as the dashed line. The spread in the predictions is due 
to uncertainties in the extraction of the GO matrix elements. A 15 % contribution 
has been added to the predictions to account for J/'if mesons originating from 
cascade decays. 
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Figure 9: Differential cross-section da/dz for 50 < PF < 180 GeV and 

Pt > ‘2 GeV. The inner error bars show the statistical uncertainty; the outer bars 
show the statistical and systematic uncertainties added in quadrature. The shaded 
band shows the prediction KZSZ (NLO, GS). The spread in the prediction is due 
to uncertainties on the charm-quark mass and on the QGD scale parameter, Aqcd- 
The solid lines show the prediction of BSW (LO, GS-hGO), where the spread in the 
prediction is due to the uncertainty on the value of the shape-function parameter. 
The dashed line shows the contribution of the GS terms only. A 15 % contribution 
has been added to the predictions to account for J/xf mesons originating from if 
cascade decays. 


30 















c 


18 


I' 16 

'6 

T3 

14 

12 

10 

8 


4 

2 


I 


■1.5 


ZEUS 


ZEUS 96-97 
KK (LO, CS+CO) 

3 X KK (LO, CS+CO) 

0.4 < z < 0.9 
Pt > 1 GeV 


-1 


t 


.... 


-0.5 


0.5 1 

rapidity (y) 


Figure 10: J/'0 differential cross-section da/dy for 50 < IF < 180 GeV, 

0.4 < z < 0.9 and pt > 1 GeV. The inner error bars show the statistical uncer¬ 
tainty; the outer bars show the statistical and systematic uncertainties added in 
quadrature. The range of the KK (LO, GS-hGO) prediction is shown as the dashed 
lines. The dotted lines show the same prediction scaled up by a factor of three. 
The spread in the predictions is due to theoretical uncertainties in the extraction of 
the GO matrix elements. A 15 % contribution has been added to the predictions to 
account for Jjif mesons originating from iff cascade decays. 
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Figure 11: J/'0 helicity parameter, a, as a function of pt for 50 <W< 180 GeV 
a) and c) 0.4 < z <1; b) and d) 0.4 < z < 0.9. The error bars correspond to the 
total experimental uncertainties. The results for the target frame are shown in a) 
and b) and the results for the helicity basis frame are shown in c) and d). In a) 
and b), the prediction of BKV (LO CS+CO) is shown as the shaded band, while 
the prediction from the BKV (LO, CS) model is shown as the dashed line. In a), 
c) and d), the data are compared with the predictions of Baranov using the GRV 
(dotted line) and KMS (solid line) unintegrated parton densities. 
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